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We study linear metric perturbations around a spherically symmetric static spacetime for general 
fiR, G) theories, where R is the Ricci scalar and Q is the Gauss-Bonnet term. We find that unless 
the determinant of the Hessian of f{R,G) is zero, even-type perturbations have a ghost for any 
multi-pole mode. In order for these theories to be plausible alternatives to General Relativity, the 
theory should satisfy the condition that the ghost is massive enough to effectively decouple from the 
other fields. We study the requirement on the form of f{R, Q) which satisfies this condition. We also 
classify the number of propagating modes both for the odd-type and the even-type perturbations 
and derive the propagation speeds for each mode. 



I. INTRODUCTION 



In cosmology, modified gravity models have been considered as a possible dynamical explanation for dark energy. 
These models have been constructed in order to give late time acceleration, however it is utterly important to see 
O ' whether such modifications of gravity will also drastically change the behavior of gravity at small scales, e.g. in our 
^ ', solar system. Among these models, some made use of the extra-dimensions [H-Q. Others implemented a similar 
$H ' mechanism by introducing a modified kinetic-term Lagrangian for scalar field non-minimally coupled with gravity 
04i3|- Most of these new theories of gravity try not to introduce gravitational ghosts in the spin-2 sector. Another 
possibility, which we consider in this paper, is to require the Lagrangian for the gravitational sector to be a function 
— I ■ of the Lovelock scalars only [lT| - [2^ , 



(N' I d^xV^fiR,g). (1) 

in 



The f{R) theories represent a subset of this general model |23l - l3l| with an interesting phenomenology |32| - |36{ . This 
. Lagrangian can be rewritten in terms of Lagrange multipliers as 

O 

A/f2 r 

S^-El d*^^[FR + ^g-U{F,0], (2) 



. where one scalar field is coupled to R, the Ricci scalar, and the other to Q = R^ — ARa/sR"^ + Ra/^fii^R"'^'^'^ , the 
' Gauss-Bonnet combination. 

: Recently some papers appeared which shed some light on this class of theories. In particular it was shown that 
on Friedmann-Lemaitre-Robertson- Walker (FLRW) backgrounds, there is a new (compared to General Relativity) 
gravitational scalar mode which propagates with a scale dependent speed of propagation. That is, its dispersion 
relation behaves like oj^ (x for large values of the wave vector k |18l - [21| . This behavior changes the effective 
gravitational constant for the matter perturbation at low redshifts. The reason for the appearance of this k^ term was 
explained in Ref. [l^l by studying another background, the Kasner spacetime. In fact, on this background, together 
with an odd mode, there are three even modes ipi. The determinant of the kinetic matrix for the even modes A 
[C 9 Aijipiijjj) vanishes in the FLRW limit. Then, in this limit, one of the fields can be integrated out from the 
Lagrangian, giving rise to a term proportional to fc**. This result shows that the number of the degrees of freedom for 
this theory depends on the background, in particular it depends on the symmetries of the background. Furthermore, 
on this same background, it was shown that, unless the Kasner manifold is very close to a FLRW one, a propagating 
ghost is always present. In fact, in the FLRW limit, it is this ghost that can be integrated out from the Lagrangian 
as its mass becomes infinitely large. 

In this paper we will discuss the behavior of the perturbation about a spherically symmetric static vacuum back- 
ground, whose metric can be written as 

ds' = glJx^dx'' = -A{r) df + ^ + +r\l- z') d^\ (3) 

Bir) 1 — 
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where z = cos9, and 9, ip are the standard spherical coordinates. So far, people have considered the way to constrain 
these theories by looking at the background solutions only (approximate or not), as this step is necessary to check if 
local gravity constraints are satisfied. However, it is important to understand, and this is the goal of this study of 
ours, whether these backgrounds are stable or not against linear perturbations, and what we can learn in terms of 
speed of propagation and ghost issues for the scalar gravitational modes. For such a theory, the background equations 
of motion read 

iBf'A' UB^f'A' WF' 2BFA' B F' A' 2F 2B F 

U = ^ 1 ^ 1 (4) 

Ar2 Ar2 r Ar A r'^ ' 

„_ 2B'^' F'B' BC 4^" FB' 2^' A' F' A' B A' FA' 6B' ^ 

R^§. (0) 

where ' stands for differentiation with respect to r. We will see that, similarly to the Kasner background, also in the 
present case in general, there is one odd mode and three even modes (under the parity transformation, 9 ^ n — 9 
and if ^ if -\- n). As for the even modes, the kinetic matrix in general has one negative eigenvalue, i.e. the theory 
does possess a ghost unless 1) the theory reduces to a subclass (to which /(i?) belongs) which satisfies the equality 
f,RB.f.gQ ~ f^ng — 0; or 2) the background manifold is more symmetrical than Schwarzschild, i.e. Minkowski or de 
Sitter. These theories then face the problem of having a propagating ghost for physical backgrounds in the scalar 
sector. It is however possible that for some models introduced to explain dark energy, this mode becomes so highly 
massive that it decouples from the relevant degrees of freedom at low energy. 



II. BRIEF REVIEW OF THE REGGE-WHEELER-ZERILLI FORMALISM 



Before studying the metric perturbation of a spherically symmetric static spacetime for f(R,Q) theories, let us 
briefly review the formalism developed by Regge, Wheeler |37j, and Zcrilli [38] to decompose the metric perturbations 
according to their transformation properties under two-dimensional rotations. Although Regge, Wheeler and Zerilli 
considered the perturbation of the Schwarzschild spacetime (namely GR), the formalism solely relies on the properties 
of spherical symmetry and can be applied to f{R, G) theories as well. 

Let us denote the metric slightly perturbed from a spherically symmetric static spacetime by — ^JJ^^ + /i^^. 
Hence h^i, represent infinitesimal quantities. Then, under two-dimensional rotations on a sphere, htt,htr and h^r 
transform as scalars, hta and hra transform as vectors and hab transforms as a tensor (a, 5 are either 9 or ip). Any 
scalar s can be decomposed into the sum of spherical harmonics as 

sit, r, 9,ip)^Yl ^^'"(^' r)Yin,{9, ¥>), (8) 

Any vector Va can be decomposed into a divergence part and a divergence-free part as follows: 

K(t,r,0,¥.) = Va$i + ^^Vfc$2, (9) 

where $1 and $2 are scalars and Eat = \/det 7 eab with jat being the two-dimensional metric on the sphere and eab 
being the totally anti-symmetric symbol with eg^ — 1. Here Va represents the covariant derivative with respect to 
the metric jab- Since is a two-component vector, it is completely specified by the quantities $1 and $2- Then we 
can apply the scalar decomposition to $1 and ^2 to decompose the vector quantity Va into spherical harmonics. 
Finally, any symmetric tensor Tab can be decomposed as 

Tabit, r, 9, ^) = VaVb^'i + -lab-^2 + \ {Ea'VcVb^i's + Eb'VcVa'^3) , (10) 

where ^i, 5*2 and vl'a are scalars. Since Tab has three independent components, ^I'l, ^'2 and ^1^3 completely specify 
Tab- Then we can again apply the scalar decomposition ([5|) to ^2 and Vfs to decompose the tensor quantity into 
spherical harmonics. We refer to the variables accompanied by Eat by odd- type variables and the others by even- type 
variables. 

What makes these decompositions useful is that in the linearized equations of motion (or equivalently, in the second 
order action) for /i^j^, odd- type perturbations and even- type ones completely decouple from each other, reflecting the 
invariance of the background spacetime under parity transformation. Therefore, we can study odd-type perturbations 
and even-type ones separately as we will do in the following. 
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III. PERTURBATION IN f{R, G) THEORIES 
A. The odd modes 



Using the Regge- Wheeler formahsm, the odd-type metric perturbations can be written as 

htt = 0, htr = 0, hrr = 0, 

hta = ^ ho,tm{t, r)Eabd''Yim{9, f ), 



hra = hi^im{t,r)Eabd''Y£m{0, (fi), 



(11) 
(12) 

(13) 
(14) 



i,7n 



Because of general covariance, not all the metric perturbations are physical in the sense that some of them can be set 
to vanish by using the gauge transformation x^^ — > +^^, where are infinitesimal. For the odd- type perturbation, 
we can consider the following gauge transformation: 



im • 



(15) 



By Aim, we can always set h2.im to vanish (Regge- Wheeler gauge). By this procedure, Aim is completely fixed and 
there is no remaining gauge degrees of freedom. 

Then, after substituting the metric into the action ^ and performing integrations by parts, we find that the action 
for the odd modes becomes 



Mil 



omitting the suffixes i and m for the fields, and 









AAi 






r 




1 


^3 


- -pi 



Aa = 



2r 



2rA'^ + 2Ai + ■' \ f - 2B' - AB^') 



2VAB 

f {f-2){AF-2BA'C) [b 
~2 



(16) 

(17) 
(18) 
(19) 

(20) 



where p — Since no time derivative of /iq appears, variating with respect to ft-o yields a constraint equation. 

However, because of the presence of /iq in the action, the constraint results into a second order ordinary differential 
equation for /iq: 



[Ai{h'a - hi)]' A3 /lo + - ^2 /il , 

which cannot be immediately solved for h^. Hence, we take the following steps to overcome this obstacle. 
Let us first rewrite the action as 



(21) 



Codd ^ Ai{hi-h'+2 



' 5 hQ + Asho - Aihi 



(22) 



so that all the terms containing hi are inside the first squared term. Using a Lagrange multiplier Q, we rewrite 
Eq. dm as follows 



^odd — Ai 



2Q[hi-h'Q + 2—] - 



'-^^^^^±^ + Mh,^^A.hi^. 



(23) 
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Then, both fields ho and hi can be integrated out by using their own equations of motion, which can be written as 

hi = — , (24) 

ho = iir A + 2Ai)Q + rAi Q'] . (25) 

These relations link the physical modes /ig and hi to the auxiliary field Q. Once Q is known also ho and hi are. After 
substituting these expressions into the Lagrangian and performing an integration by parts for the term proportional 
to Q' Q, one finds the Lagrangian in the canonical form 

where 

Air"^ { r'^A\ A' - r^A'.'As + 2 A1A3 + 4A\^+ As^r^ - 2 AiA'' + 2 Air A' - 4 A\ rA^ 

u2 ^ 1^ t : t t : ^ (27) 

{2 Al + 2 A[r - Asr-^ f 

From Eq. (j26p . we can derive the no ghost condition 

^4 > , or equivalently AF -2B^' A' >0. 

For solutions proportional to e*^"*"'''') with large k and uj, we have the radial dispersion relation 

2 _ B {AF -2BC A') 2 
~ (F - 4Be" - 2B'e') ' 

where we made use of the background equations of motion. Finally the expression for the radial speed reads 



-"odd 



dnV {AF ~2B£,' A') 



dr J A{F- 45^" - 2B'^') ' 



where we used the radial tortoise coordinate (c?r^ = dr^ / B) and the proper time {dr'^ — Adt^). Therefore in order 
for the modes to be stable, one also requires 

F - ABC - 2-8'^' > . 



B. Even modes 



Now that we have got an idea how the action approach works for the odd modes we can tackle the more complicated 
problem of the even modes. In this case, the perturbed metric can be written as 

htt = -A{r) J2 HoMt^ r)Ye^{9, <f), (28) 
htr = Hi,tm{t, r)YiUO, (29) 

£,171 

hrr = -Jj--y^H2Jmit,r)YijniO,f), (30) 

hra ^'Yaemit,r)daYi.miO,(p) , (31) 

and we use the gauge transformation to set hta and hat to vanish. In addition to the metric perturbations, we need 
to perturb also the extra scalar fields F and ^ as 

F = Fir) + J2SFernit,r)Yim, and ^ = + 1] r)y,„, . (32) 
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Then, the action at second order for the even modes, reads as follows 

Seven — ^ ^ ^ / dv /Icvcn 5 (^^) 

where 

£even = {aiS^' + a2 SF' + a^H2 + a4 a + a5 (5^ + 5F) + (a? + as a + ag + aio 

+j2 61 iJi^ + (62 5^ + 63 +biH2+ f b5a + beS^ + bjSF) + ci i?| 
+i?2 [C2 SC + C3 (5^^' + f c^a + S^ if C6 + cj) + SF [f cg + cio)] + (C5 k + C8 5>) 

(di + d2 a^) + f a (dg 5^' + ^4 + ^5(5^ + 4 5F) + ei + 62 Ji^ 5^ + /i , (34) 

where ai,bi,Ci,di, Ci and /i are all functions of r only and their expressions are given in the Appendix. For simplicity, 
we omitted the subscripts m also for the even modes. In what follows, we will integrate out the fields Ho, Hi, and 

We first integrate out the non-propagating field Hi , by using its own equation of motion 

Hi=- -^T- {b2 k' + h SF +b4H2+fb5a + be 5^ + bj SF) . (35) 

We note that the term proportional to Hq is absent in the action. Thus, the equation of motion for Hq 

ai S£," + a2 6F" + (05 + a[) S£_ + (4 + ag) SF + f a' + as H^ + - f ag) S^ + {a'^ - f aw) S^ 

+ (4 - f ar) H2 + « - as) « = , (36) 

sets a constraint for the other fields. By looking at this equation, one may think that we cannot use it to directly 
substitute back any of the fields. However, we can remove the highest r-derivatives for the fields SF, S^, a and H2, 
by performing the following field redefinition 

a/ia — a^vo — 03 H2 — a2v'i . (37) 

SF = vi-^S^^vi-^^^^SC (38) 
02 

Because of this field redefinition, Eq. p6p has no more second r-derivative for any of the fields, and no r-derivatives 
for the field H2. Solving for H2, we obtain 

3 aj -^2 = ae H Vi + I ar^ - a^ ^ dS, +a4V„ + (og - oiq) vi + (04 - ag) Vq 

0,4 J V a4 / \ a2 a2 I 



2 „ , «i (4 - ag) , ai (a^' - a^) ^ ai a'2 (og - a^) ^ p ai a^ 



02 a2 a, a2 



SC. (39) 



We can now substitute H2 back into the action, so that also the field Hq is automatically eliminated. We will find it 
convenient to finally perform the field substitution 

«o = t'2(l + 4j2)i/2^ (40) 

SC = v,il + 4f)'^', (41) 

where the dependence on j in Eqs. (|40p to (1411) is chosen such that, for large j, the ghost conditions, which will be 
found below, become independent of j. Now the Lagrangian takes the canonical form 

3 

Ccvon = ^ [Kij ir, j) Vi Vj - Lij (r, j) v[ v'j - Dij ir, j) v[ Vj - Mij (r, j) Vi Vj] , (42) 

ij = l 

where i,j run from 1 to 3, and the coefficients important for our discussion are given in appendix. All matrices are 
symmetric except for Dij, which is anti-symmetric. Now we can discuss the existence of ghosts and the speed of 
propagation for the modes. The no-ghost condition requires the matrix K to be positive definite, that is 

K33>0, K22 K33 ~ Ki:, > , det(i^,,)>0. (43) 
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We find that K22 -f^33 — -f^ls given by 

;(l + 4j2)2AB(2B- 

^2 ^2 



16 (1 + Aff AB{2B-2-r B'f [4(i? r^Ff ^ ^ 

iV22A33 - i<23 = j-r^ ^U, (44) 



where 

A = 24AB2^' - UB^A'^' - 8AB^' + Ar^'BA' - AABFr + 2BFA'r^ - 2ABF'r^ + BA'r^F' + f {2AFr - SAB^) ■ (45) 
Therefore, on this background, a ghost is always present. The determinant of the kinetic matrix is given by 



if - 2)(1 + 4j2)2j.2(2S -2- rB')^ {F - 2B'e - 45^") 
-2^2 



det(if.,) = ~32VAB '-^ -A^^^7 , ' ^ i ' . (46) 



This quantity gives us new information. Indeed, on backgrounds of exact solutions given by i? = 1 + C with C being 
constant (Minkowski or de Sitter solution), both the determinant and the ghost kinetic term vanish, which implies an 
effective reduction of the degrees of freedom similar to what occurs on FLRW background, where the missing degree 
of freedom is ghost-like [l^. The results obtained for the odd modes impose F — 2i?'^' — ABSf' > 0. Then, for the 
even-modes, det{Kij) < holds, so that either there is one ghost (X33 > 0) or there are three ghosts (-fi'33 < 0) on 
the backgrounds of general spherically symmetric static solutions in these modified gravity theories. The cases with 
£ = and £ = 1 should be treated with care, as some variables are absent from the beginning due to the lack of vector 
(for ^ = 0) and tensor (for ^ < 1) harmonics: for example, there is no contribution from a in the action for j = 0. 

Assuming the solutions in the form proportional to e"^*"''^'", the dispersion relation for large k and w can be obtained 
by solving the discriminant 

det(cj2 Aij ~ fc2 Dij) = , (47) 

which is a cubic equation in w^. The three radial speeds of propagation we obtain are 

^ ^ {2AB-2A~rBA') 

i2B-2-rB')A ' ^''^ 
„2 AF-2A'Be (^g) 



A{F-2B'^' ~AB f ) 



Two of the speeds of propagation reduce to unity for backgrounds with A = B. The third one, which is identical to 
the one for the odd modes, depends directly on the profile of the two new scalar degrees of freedom F and ^ even if 
we set A = B. The speed of propagation for large j, does not have a simple analytical form. Although one can solve 
a cubic equation in w^, the expression is too complicated to gain intuition from it. 



C. Discussion regarding the ghost 

We have found that at least one ghost mode is always present around the spherically symmetric static background 
in vacuum. However, the existence of the ghost does not necessarily mean that the background spacetime is unstable 
due to the creation of the ghost and normal particle pairs. Our implicit assumptions are that there is an yet unknown 
complete theory (maybe string theory), which is well-defined at any energy scale and does not have any ghost and that 
the f{R,G) theories are the derived effective theories, which are valid only below some cutoff scale Mcutoff- What is 
generally thought is that if the mass of the ghost is always heavier than Mcutoff, such a ghost should not be regarded 
as a physical mode and must be integrated out to have more sensible effective theories. On the other hand, if the 
mass of the ghost becomes lighter than Mcutoff in certain situations, such theories do not make sense and must be 
ruled out from the list of the possible low energy effective gravitational theories. 

In the models we study here, there are in general two ways out of which the ghost mode can become massive. 
One is due to symmetry. That is, as the background becomes more and more similar to either Minkowski, de Sitter 
or FLRW, its mass tends to infinity because its kinetic term vanishes. The other is due to the so-called chameleon 
mechanism. Thanks to the local value of R or G much larger than the corresponding cosmological value, some modes 
may develop their masses. Our cutoff mass Mcutoff is not necessarily as large as Planck scale, but M^^^^^^ must be 
sufficiently smaller than the experimentally relevant length scale Lcxp- 

Let us give approximate values for the masses of the modes, including the ghost mode. Assuming the background 
around a star to be very close to the standard GR case, then one has F ss 1, or F' w 0, and ^' w « Further we 
assume that the theories satisfy solar system constraints, that is, A w i? sa 1 — r^/r, (where rg is the Schwarzschild 
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radius of the star). Under these assumptions, the leading contribution in the mass matrix is given by the terms 
originating from U,ff, U^fs,, and f/.j^, 



1 9 



(1 + 4j2)i/2 u^pp _ J.3 u^p^ 



M: 



1 + Aj' 



33 



2 r 
Then the discriminant equation, for low fc, to solve is 

det(TO2 A{r) K^j - M^j) = , 



2r4 



FF) ■ 



(50) 



where the factor A{r) comes because of the choice of the proper time as time variable, and the elements Kij are given 
in Appendix. Equation (I50p reduces to 



gg + 3(C/.5^ - GU^ff) to^ - {U^FF - U'^FV 



= 0. 



(51) 



where we introduced Q = 4t?/3 ~ 16r^/r^. It is also suggestive to rewrite the above equation in terms of the function 
/(_R, Q), by using the relation 



as 



m 



U^FF U^Fi ] _ ( f,RR f.RS 
U^FS, C^,?? / V ^'R^ ^^^^ 



90 det(/,,,) + 3(/^^^^ _ gf^^^) ^2 _ ^ 



0. 



(52) 



(53) 



One obvious solution of Eq. ([ST|) is = 0, whose eigenvector is given by Eq = [0,1,0]*. The kinetic term of this 
mode 



E^AEo 



P [3r, + i f - 2) r]2 



is positive for £ > 2. Thus, this mode represents the GR standard contribution. The other two mass eigenvalues are 
given by 



^U,, + GUff ± {U,, + QUff)^1- j^l,^^,. 
i| = ^ 

eg 



(54) 



In order to study which of these mass eigenvalues correspond to the ghost mode, let us perform a little more detailed 
analysis. First, we diagonalize the kinetic matrix A Kij by arranging linear combinations of fields Vi = Pij Wj, which 
are explicitly expressed as 



vi = 



V2 



V3 



Vq 2 



W2 , 



^ / -wi, ^ , 

■ W2 ■ 



r j + 7rs — 6r 



W2 +W3 . 



2V6 yrrifi 

Then, the new diagonalized kinetic matrix Kij = A{r) Pki Kki Pij takes its diagonal elements 

2(r-r,)2(l + 4j2)(j2_2) 



(55) 
(56) 
(57) 



i^ii = 1 , K 



22 



-1, K 



33 



f [{f -2)r + 3r, 



whereas the new symmetric mass matrix, Mij = PkiMkiPij, satisfies M^s — for i = 1,2,3, and its non-zero 
components are 



U 



FF 



U 



Mil = , - _ 

3 sVg 



M22 - 



u 



35 
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At this point, one can easily see that M22 is the most enhanced mass-element in the Minkowski limit, — 0, which 
gives a divergent mass to the ghost mode. This is consistent with the notion that the ghost possesses a divergent 
mass in the Minkowski background. However, one possible natural hierarchy among the components of Mij will be 
that Mil, and M22 are the same order when Q = 0{Hq), where Hq is the present value of the Hubble parameter. 
In this case, since the value of G around the local gravitational source is typically much larger than the cosmological 
value Hq, the simple limit G 0, which relatively enhances M22, is out of the relevant parameter range. 

Since the GR mode completely decouples from the others, we shall concentrate on wi and W2 below. They 
are still coupled through the off-diagonal mass matrix Afi2. The field transformation which keeps the kinetic matrix 
unchanged is given by wt — Zijzj, where Zn = Z22 ~ cosh/3, and Z\2 = Z21 = sinh/3. Then, the condition for the 
off-diagonal component of the new mass matrix to vanish yields 

tanh2/3 = -^^^?^. (58) 
U^^ + GUff 

Therefore, AGU^p^ < {U^^^ + GU,ff)'^ is required to have the mass matrix diagonalized. In terms of /, this condition is 
equivalent to —AG det(/_y) < {f^Hri. — Gf,gg)^- In this case, the new fields can be identified as independent decoupled 
massive modes, and 11 and 22 components of the diagonal mass matrix are identified with and —rn?_, respectively. 
Namely, and m?_ are the squared masses of the non-ghost and ghost modes, respectively. This identification can 
be easily verified by considering the trivial case with t/,F{ = 0. The condition for the model to be applicable to an 
experiment on a scale Lexp will be |m^| ^ -M^^^^jj, with M^^^^g > L~^^. If we do not see any significant deviation 
from general relativity on this scale, a condition |m^| > has to be imposed as well. These conditions are less 
intuitive. If we are allowed to crudely identify M^^^^jj with , the above conditions are simplified to 



9^det(/,, 



Hf,B.R - Gf, gg) 





1 


1 










771 q_ 





« icxp , (59) 



where we used the fact that 771^ and rri^ are solutions of Eq. (|53p . 

Next we consider the case with 4^[/^^ > {U^^^ + GU^ff)'^ (or equivalently — 4^det(/,.y ) > (/_i{fl — ^/_gg)^), in which 
we cannot simultaneously diagonalize both the kinetic and mass matrices. In this case the eigenvalues m? obtained in 
Eq. ([54]) become complex. The complex nature of the solution means that those modes are classically unstable. Since 
the eigenvalues for 771^ are complex conjugate with each other, a unique mass scale y^|r77?| must be much larger than the 
cutoff scale. Again, with the aid of Eq. ([5^ . this requirement leads to -~{G det/^j) < L^xp- This condition combined 
with — 45 det(/.ij) > {f.RR—Gf,gg)^ is identical to the conditions of the previous case, given in (|59|) . Thus, we conclude 
that the conditions ([5^ are the necessary and sufficient conditions for the model to be viable. When we apply these 
constraints on the model to the solar system, we need to plug in i? « (47r/3)yOiocai ~ (47i'/3) x 10^^'*g/cm'^ « (lO^^cm)^^ 
and G ~ 48MQ/r^ ^ (5.7 x 10"'^^cm)~'''(lAU/r)^ as the background values. In the following we present two examples, 
in which the expressions are simplified by assuming some hierarchy among the components of the mass matrix Mij. 

The first case is the one in which the off-diagonal element M12 is suppressed, i.e. f^ag is negligible. If we set 
Afi2 — 0, the expressions for the mass eigenvalues simplify to give 

2 ^ U,FF _ 1 2 _ _ 1 frn\ 
771 I « — - — « — , m_ « 4^ « = . (60 

+ 3 3/,flfl' 3^ 3Glgg 

The conditions that these masses are positive, f^B.R > and f gg < 0, imply det(/ jj) < 0. Such a situation is realized 
by considering the following type of toy models of dark energy: 

fiR,G) = R + TRiR)+Tg{G), (61) 

where the functions J^fl(_R), and J-g{G) have to be chosen such that /,_r_r > and /,gg < 0. The correction J-r, is of 
the kind discussed in Refs. [28l. l29l. Isol. l40l|. whereas Fg is similar to the functions introduced in Refs. [13, El. For 
concreteness, we specify the functions as 

FR{R)=AR^l''P+V{RP + CR^i^'n , -Fg(e) =-4g/"+V(^2„_^^^^8n) ^ 

where /i = 0{Hq) and AR,Ag,CR,cg are constant parameters of 0(1). In these models, the corrections to the 
cosmological evolution become important only at around the present epoclQ. When we consider local gravity in a 



^ Here, Tg is a good approximation to the models introduced in |14|| when Q ^ H^, as it happens in the vacuum Schwarzschild solution. 
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dense region, the values of R and G are much larger than the cosmological backgrounds: R ^ Hq and G ^ Hq. In 
this case, we have U^ff = = 0{fi-^P~^ RP+^) and U^^^ = l/f,gg = 0(m"*"-2 g2n+2-)^ ^^^^^ ^j^g j^^^^g squared 

for each mode is given by 

'"■non-ghost ^ ^ I ^0 I J '"-ghost 



Pc 



o m 



2 ,4n+2N 



(62) 



where ~ 4 x 10~^'^g/cm'^ is the critical density of the universe. For the non-ghost mass, we have used the Einstein 
equation, which is a good approximation in the local region like the solar-system, to replace the Ricci scalar with 
the energy density piocai of matter surrounding a staiQ. If the system is exactly the vacuum, the non-ghost mode 
becomes massless. Therefore, for the toy model of Eq. (|61]) . the non-ghost mode acquires mass through the chameleon 
mechanism. For example, if we put the values for the solar system, Ts ^ 3km, r ~ lAU, piocai ~ 10~^^g/cm'^, then we 
have '77,j^jQjj_gjjQg^ « (10^^ ''^^ ''Pcm)^^ and m^Yiost ~ (10^~®* ''"'cm)^^. One can easily make the ghost mode sufficiently 
massive, while we need to choose a relatively large power p > 4 to make also the non-ghost mode sufficiently massive, 
i-e- <o„-ghost » (lAU)-2. 

Interestingly, unlike the standard chameleon mechanism, we do not need the matter to make the ghost mode very 
massive. However, the background value of G is not always much higher than the cosmological value. For example 
the value of the Gauss-Bonnet term inside a star is typically negative (if evaluated, e.g. on a Schwarzschild interior 
solution) whereas it will be positive outside the (neutron) star. Therefore there should be a point where G switches 
sign. At this point, the light ghost problem might arise. However, since the region of the surface of the star where 
G = 0{ii^) might be very thin, the notion of the mass in the current treatment, in which spatially homogeneous 
modes are assumed, becomes irrelevant. Hence, it is not clear if actually there is a problem of light ghost at this 
point. Another way around this problem consists of thinking of a theory in which the background value of R also 
contributes to giving mass to the ghost when G is small. Some modification of Fg like J-g cx l/(tj^" + ci?'' + cg/i*") 
with a large power q may work, but still some engineering for choosing the function will be needed. 

The second case is when \G det(/,y )| ^ {f,RR — Gf,ggY- Basically, this is the case when the matrix f^ij is almost 
degenerate. Once the condition \G det(/^ij)| ^ {f,RR ~ Gf, gg)^ holds, the mass eigenvalues reduce to 

"^non-ghost ~ ^J— : Wghost ~ " „^ , ... . ' for \I M > l^/.Sel, 

^{f,RR- Gf ^gg) 3^^det(/,y) 



f,RR - Gf, gg 2 ^ 

3^det(/,y) ' ^ 3(/,i^fl - g/,ge) 



"i-non-ghost ~ ^ ''oXj ^/j 'T ' '"ghost ~ —77 T 7 for \Gf, gg\ > 



and the mass hierarchy iTOgj^^^J ^ I'^non-ghostl automatically guaranteed for > \G.f. gg\- In contrast, when 

l^/.eel > 1/,-R-rI; the mass scale of the ghost mode is always lower than that of the non-ghost mode. This limit 
includes a natural situation in which f,RR, iJ^f.Rg and fi^f.gg are the same order with the typical energy scale fi 
satisfying ^ G- In this case, we have \Gf, gg\ ^ |/,fli?J and therefore the mass scale of the ghost mode is lower 
than that of the non-ghost mode. A more explicit expression for the masses are 'rn^^^_^Y^^^^ ~ det(/_ij)) and 

"^ghost ^ ~1/(3^/,Gg)- The expression for 'rn^^^_^^^^^ is slightly different from the previous example in which we 
simply neglected f^Rg- 

D. Case £ = 0. 

In the case with j = neither tensor nor vector harmonics exist. Thus, the field a does not contribute any longer 
to the action. Furthermore also the term quadratic in the field Hi disappears, together with the linear term in Hq. In 
this case the variation with respect to Hi gives a constraint which can be solved for H2 in terms of SF, 5^, 6F' and 
6C- At the same time also the contribution from Hq in action (p4)) automatically cancels. Therefore we are left with 
an action in terms of SF, (5^ only, i.e. there are only two degrees of freedom. In this case we find that the determinant 
of the new kinetic matrix A^'^^^ can be written as follows 



det 



\ '3 ) AB (12Be - 2rF - 4^ - r^F')^ ' ^ ' 



2 



Strictly speaking, one cannot take into account the effect of the matter piocal in the present discussion since we have assumed vacuum 
throughout the analysis. However, if the gravity is dominated by the central star, which is true for the solar-system, we expect that the 
gravity from the surrounding matter can be neglected except for shifting the background values of R in evaluating the matrix elements 
of U^ij or f^ij, and our vacuum results can be used as the first approximation. 
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Therefore also in this case there is one (and only one) ghost degree of freedom. The speeds of propagation for both 
modes are equal to the expressions given in Eq. 



E. Case £ = 1. 

In the case with = 2 the determinant in Eq. (j46|) identically vanishes, as the tensor harmonics do not exist. This 

(£—1) 

leads to a reduction of the propagating degrees of freedom. In fact now the matrix A]- possesses an eigenvector 
N2 with zero eigenvalue which satisfies A^^^-'^^ N2 = 0. This vector can be written as 

3(45^' - 4^' - r^F') " ^ * 



[ r2f ' 2AB(rF-4Se)'7 ' ^^^^ 

where 

T = -4:ABFr + 2BFA'r^+BA'r^F' + 4:FrA 

- 12B^^' A' r + 32AB^^' -8^' BArB' -S2AB^' + 4:A' Br^' . (65) 

This suggests the field redefinition 

3(45^' - 4^' - r^F'^ 
72^ 



vi=Qi YT/ ^' (66) 



"^=^^+2l5(7y3i^^' (67) 
V3=X. (68) 
In this way also the couplings between X with Qi, and vanish. The Lagrangian reduces to 

C = Q, - bff^^ Q[ Q; - M^'=') Q, + Ci + C,X , (69) 

(with i, J — 1, 2) so that by integrating out the field X, only the mass term is affected. In this case we can find the 
ghost condition and the speed of propagation. In particular 

and therefore also in this case one of the two propagating degrees of freedom is a ghost. The speeds of propagation 
for the two modes are equal and coincide with the expression given in Eq. (j^S]). 

IV. SPECIAL CASES 

In the theories of gravity discussed so far, we have studied the general function /(-R, Q) for which 

0. (71) 



dR^ ag^ \dRdg) ^ 



However there is a special class for which s identically vanishes, e.g. the f{R) theories. In general, S vanishes when 
F and ^ are not independent, i.e. when ^ = £,{F). When this happens, we have 

S^^jjSF, (72) 

and the independent scalar degrees of freedom reduce by one. In this case one has to follow the same procedure as 
we have done for the case S 7^ 0. In particular now 

= ^'i , (73) 

and a is substituted by vq according to 

J a4 a = a4 Wo - 03 ii2 , , „, Wi 

Abi }'' 
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by using the same relation (j37p . In order to remove higher r-derivatives from the action one performs another field 
redefinition as 

^ = «i , (74) 
vo = il+ 4j") V2 . (75) 

We have two no-ghost conditions which must be satisfied in order to remove any ghost degree of freedom. In particular, 
the determinant of the 2x2 kinetic matrix Aij becomes 

if - 2) (1 + Af)^ 4A B [r' F' - A{B ~ 1) (F - 2B' C - 45^) T2 
det(^u) = -2 {F')^A^ ' ^ ^ 

where 

r2 = 3 r^FF' + 2 ^'r^F'B' - 4 ^'F - 16 C'B^C + 4 Br^F'^" 

+ 16 ^B^" + 24 C'^BB' - 8 r^'FB' + 8 ^'^B' - 16 C'BF'r + 4 ^'BF , (77) 

and A is defined in Eq. (|45|) . The other independent condition, say An > 0, is rather complicated, but for large j is 
given by 

, 2irF' -2B'C')[r^F' -4{B-l)e] „ 

lim All = — , ' ^ — > , (78 

j^^ ^/AB {F'Y {r F - AB i') ~ 

so that in this case the ghost can be absent. As for the speeds of the two modes, one coincides with the expression 
given in Eq. P5|) . whereas the other becomes 

where 

Eg = 24 C'^B^ A' + 8 C'^BA' -Ai'F A + 2 ^'r^F'A'B - 8 ^'rA'F B - 16 C'F'rAB + 4 ^'FAB + 3 F'r^FA . (80) 



V. CONCLUSION 



We have studied linear perturbations around the static spherically symmetric spacetime for the modified gravity 
theories whose Lagrangian consists of a general function f{R,Q) of both the Ricci scalar R, and the Gauss-Bonnet 
scalar, Eor the odd- type modes, there are two degrees of freedom (one dynamical variable). We derived the no 
gradient instability condition. These conditions put constraints on the background quantities. 

For the even-type modes, the picture is more interesting. We have found that there are, in total, four degrees of 
freedom (corresponding to two dynamical fields) for monopole and dipole perturbations {£ = 0, 1), which do not exist 
in GR. Both for the monopole and dipole perturbations, one of these two new scalar modes is always a ghost. As 
for the higher multipole perturbations, there are, in total, six degrees of freedom (corresponding to three dynamical 
fields), four of which do not exist in GR. Eor Minkowski and (anti-)de Sitter backgrounds, one out of these three 
fields is not independent from the others. Instead, for general backgrounds (including Schwarzschild) , a ghost always 
appears from the even- type perturbations with mass m_ given in Eq. (I54[) . whereas m+ is the mass of the non-ghost 
chameleon field. Einally the third dynamical field is massless describing the standard GR-mode. The necessary 
and sufficient conditions that both the ghost and chameleon fields are sufficiently massive, i.e. their inverse mass 
scale is smaller than iexp, the length scale of a determinate experiment, reduces to ([5^ . When these conditions 
hold, the ghost is massive enough to be treated as a non-propagating mode in the effective gravitational theory. In 
order to satisfy both these conditions, some hierarchy among the fjj coefficients is required, which means that some 
engineering in choosing the function f{R,g) is needed. 

We also found that, in the high frequency limit, the radial propagation speed of one field, among the three even- 
type modes, coincides with the one for the odd-type mode, and the remaining two even-type modes have a common 
propagation speed which is different from the former. The classification of the modes are summarized in Table I. 



^ These theories are equivalent to the double-scalar-tensor theories defined hy C = FR + ^C} — U{F,^). 
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odd-type G 


odd-type S 


even-type G 


even-type S 


number of modes 


(for/ = 0, 1) 
1 (for / > 2) 


(for £^0, 1) 
1 (for £>2) 


2 (for £ = 0, 1) 
3 (for £ > 2) 


1 (for £ = 0, 1) 
2 (for £>2) 


ghost and 
gradient conditions 


AF - 2A'BC > 
F - - 2B'^' > 


AF - 2A'B^' > 
F - 45^" - 2B'^' > 


ghost present 
massive if Eq. (I59II is verified 


constraints 
(see SecHVt 



Table I. Classification of the modes for the general f{R, Q) theories. G and S stand for the model that does or does not satisfy 
Eq. ISIl). 



If the theory satisfies the condition: 



dR^dG^ \dRdG ' ^ ' 



then the number of modes for the even-type perturbations reduces by one. In this case, the ghost is not necessarily 
present, i.e. we get non-trivial no-ghost condition which puts bound on the function form of f{R, G) and also the 
background spacetime. For example, in /(i?) gravity theories, the no-ghost condition is satisfied once df /dR ~ F > 0, 
and we find that the ghost does not exist in this case as expected. 
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Appendix A: Coefficients in the action 



We define here the different coefficients introduced in order to define the action. By calling 6 = \J AjB, we have 



ai = 


-04 t)2/(46i) , 


(Al) 


OjI = 


n . /lo zL/it ^ 

"4 "3/ J ) 




rto 

"3 — 


"4 "4/ v^i'l y ) 


(A3) 




B R Cr F — 4 /r^ 


(A4) 


05 — 




[AO) 








Oil = 


"-^^t '-'0 / "^.^ 5 


(A7) 


as = 


-65 (4 B^' - F r) (A' r - 2 ^) /(2^r^) , 


(A8) 


(Xg = 


-2eB7r^, 


(A9) 


dlO = 


Cg , 


(AlO) 


61 = 


QB (Fr-4BC) /(2Ar^), 


(All) 


62 = 


8Be{B-l) /(Ar^), 


(A12) 




2QB 


(A13) 


03 = 


A ' 


64 = 


-QB {-F'r^ -4^' -2Fr + 12B^') /{Ar^), 


(A14) 


&5 = 


-2 61, 


(A15) 


66 = 


&2 br/bs , 


(A16) 


67 = 


-46106/04, 


(A17) 


Cl = 


-Se (24 B A' - 4 ^'A' - r^i^^' ^' - 2 F r A' - 2 ^ - 4 ^ r F') / (4Ar2) , 


(A18) 


C2 = 


2eBA'{3B-l)/{Ar'^), 


(A19) 


cs = 


-Be {rA' +4:A) /{2Ar) 


(A20) 


C4 = 


BQ (12 B^'A' -2ArF' -2AF - FrA') /{2Ar^) , 


(A21) 


C5 = 


4e(S-l)/(^r2), 


(A22) 


C6 = 


da C5/62 , 


(A23) 


C7 = 


-e{BA''^-2ABA"-B'^A'^ + 3ABB'A' + 2B^AA" -AB'A')/{A'^r'^), 


(A24) 


C8 = 


&3C5/&2 , 


(A25) 


cg = 




(A26) 


ClO = 


e{AA^B' - BA''^r + 2ABA"r + AB'A'r) / [ArA^ ) , 


(A27) 


di = 


bi, 


(A28) 


d2 = 


-QB {2B£_'A'-A F) /(Ar^) , 


(A29) 


d3 = 


-4eB^A'/{Ar^) 


(A30) 


d4 = 


Cg 62/C5 , 


(A31) 


dr, = 


AGB^A'/iAr"^) , 


(A32) 


de 


-2Se/r^ 


(A33) 


• = 


-ef/,FF/2, 


(A34) 


62 = 




(A35) 


/l = 




(A36) 
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The elements of the kinetic matrix relevant for our discussion are 

K22 = -tTa r~2~^ nm ['^^l^l iha'sai + -2bial - b'^asaA 

— Sj^bltti (a4C562^'4 + 040562^4 — a4&2C5&4 — 2 0562^404 + 2 C5a8''2&4)] (A37) 

8j2 6f 6|r (4050461^2 + a86264)2 L ' ^ ^ ^ 

— 2 j^6f (2 a\czhirh'2 — 4 a/^c^b2bia'^r + 2 alc5b2b'^r — 62^64^0809^ + 8 a4C5a86264r — 2 a|6264C5r) 

+ 6|64r (6ia8a4 — 2 6103 + 610403 — 040861)] , (A38) 



where 



b2b'3-b3b'2(xrB' -2B + 2, 



which vanishes on (anti-)de Sitter and on Minkowski. 

For backgrounds close to the Schwarzschild one, the matrix Kij reduces to 



ifii 
K12 

^13 



(2 r'^f + (10 rs - 7r) + (14 r^^ - 20 rr^ + 6 r^) f - r^^) 



j2 ^j.j2 _ 2 r + 3 fg)^ (r — Tg) 



y.2 



v/rT47(j2-2) (rj^+r,) 



j2 (^j,j2 _ 2 r + 3 Ts) 
6 (r^i*^ + 2 (2 Ts - r) rj"* + (7 r^^ - 6 rr^) - 4 rr^ + 4 r^^) yTTT^ 
j2 (-j,j2 _ 2 r + 3 Ts)^ (r - rs) 
2(r-r,)r(l + 4j2) (j2 _ 3) 



iMri2-2r + 3r«)' 
-K'23 = K22 , 

-f^33 = — TT- K23 ■ 
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